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Abstract 

Background: NAC domain transcription factors are important transcriptional regulators involved in plant growth, 
development and stress responses. Recent studies have revealed several classes of NAC transcriptional factors 
crucial for controlling secondary cell wall biosynthesis. These transcriptional factors mainly include three classes, 
SND, NST and VND. Despite progress, most current analysis is carried out in the model plant Arabidopsis. Moreover, 
many downstream genes regulated by these transcriptional factors are still not clear. 

Methods: In order to identify the key homologue genes across species and discover the network controlling cell 
wall biosynthesis, we carried out comparative genome analysis of NST, VND and SND genes across 19 higher plant 
species along with computational modelling of genes regulated or co-regulated with these transcriptional factors. 

Results: The comparative genome analysis revealed that evolutionarily the secondary-wall-associated NAC domain 
transcription factors first appeared in Selaginella moellendorffii. In fact, among the three groups, only VND genes 
appeared in S. moellendorffii, which is evolutionarily earlier than the other two groups. The Arabidopsis and rice 
gene expression analysis showed specific patterns of the secondary cell wall-associated NAC genes (SND, NST and 
VND). Most of them were preferentially expressed in the stem, especially the second internodes. Furthermore, 
comprehensive co-regulatory network analysis revealed that the SND and MYB genes were co-regulated, which 
indicated the coordinative function of these transcriptional factors in modulating cell wall biosynthesis. In addition, 
the co-regulatory network analysis revealed many novel genes and pathways that could be involved in cell wall 
biosynthesis and its regulation. The gene ontology analysis also indicated that processes like carbohydrate 
synthesis, transport and stress response, are coordinately regulated toward cell wall biosynthesis. 

Conclusions: Overall, we provided a new insight into the evolution and the gene regulatory network of a 
subgroup of the NAC gene family controlling cell wall composition through bioinformatics data mining and bench 
validation. Our work might benefit to elucidate the possible molecular mechanism underlying the regulation 
network of secondary cell wall biosynthesis. 
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Background 

As a potential replacement for traditional fossil fuels, bio- 
fuels have received increased public and scientific atten- 
tion in recent years [1], The current first generation 
biofuel is based on sugar and starch derived from feed- 
stocks such as sugrarcane and corn; however, this platform 
is not sustainable for various reasons. Lignocellulosic bio- 
mass has been proposed as the major feedstock for second 
generation biofuels to enable the transition from fossil 
fuel-based energy to renewable energy for the various eco- 
nomic and environmental advantages gained over first 
generation biofuels [1,2]. Generally speaking, lignocellulo- 
sic biomass is composed of cellulose, hemicellulose, pectin 
and/or lignin, but the amount and ratios between the 
components can vary considerably [3]. In addition to the 
aforementioned components, even the amorphous por- 
tions of cellulose are purportedly important for lignocellu- 
losic conversion to biofuel [4-7]. Since secondary cell walls 
in fibres and tracheary elements constitute the most abun- 
dant biomass produced by plants, it is necessary to eluci- 
date the possible molecular mechanisms underlying the 
regulation of secondary cell wall biosynthesis for improved 
plant biomass production. 

Plant NAC (NAM, ATAF1/2 and CUC2) domain pro- 
teins are one of the largest groups of plant-specific tran- 
scriptional factors and are known to play diverse roles in 
various plant development processes and stress response. 
NAM (no apical meristem) was the first characterized 
NAC gene in petunia. The NAM gene product is 
required for apical meristem formation and correct posi- 
tioning of the cotyledons during petunia embryogenesis 
[8]. ATAF1 and ATAF2 are the two NAC genes in Arabi- 
dopsis playing negative roles in response to drought and 
pathogen infection respectively [9,10]. CUC2 (CUP- 
SHAPED COTYLEDON 2) gene was also characterized 
as a NAC gene in Arabidopsis [11]. Arabidopsis RD26 
(RESPONSIVE TO DEHYDRATION 26) encodes a NAC 
domain protein [12] with function in ABA-mediating 
gene expression under stress conditions [13]. StNAC, 
one potato NAC gene, was shown to be rapidly and 
strongly induced by wounding [14]. Over-expression of 
OsNAC6/SNAC2 in rice can enhance the seedling plants 
tolerance to drought, salt, and cold stresses [15,16]. 
Recently, accumulating evidence has indicated that a 
considerable portion of NAC domain proteins play cru- 
cial roles in the processes of xylogenesis, fibre develop- 
ment and wood formation in vascular plants [17]. In the 
model plant Arabidopsis, NST1 (NAC Secondary Wall 
Thickening Promoting Factorl), NST2 and NST3/SND1 
(Secondary Wall-associated NAC Domain Proteinl) are 
key switches in regulating secondary cell wall biosynth- 
esis in a partially redundant manner [18-25]. NST1 and 
NST2 function redundantly in regulating secondary cell 
wall thickening in the endothecium of anthers whereas 



NST1 and NST3/SND1 were shown to regulate second- 
ary wall thickening in fibres. In Medicago sativa, 
MtNSTl (Medicago truncatula NAC Secondary Wall 
Thickening Promoting Factor 1) has been identified as 
the only homologue of AtNSTs [26]. Loss of function of 
the single MtNSTl gene resulted in lack of lignifications 
in interfascicular fibres, loss of anther dehiscence and 
stomatal phenotypes associated with loss of ferulic acid 
in guard cell walls. VND6 (Vascular-related NAC- 
domain6) and VND7 are key regulators in protoxylem 
and metaxylem development. VND6 is specifically 
expressed in the metaxylem of Arabidopsis primary roots 
whereas VND7 is expressed in the protoxylem. Recently, 
the VND6 gene was discovered to regulate xylem forma- 
tion by directly targeting some genes related to secondary 
cell wall formation. VND6 also acts as a direct regulator 
of genes related to programmed cell death [27]. SND2 
and SND3 were also found to function in the formation 
of secondary cell walls in fibres, and were down-stream 
of NST1 and NST3/VND1. Six NAC genes associated 
with wood formation in Populus were also reported [28]. 
Among the six genes, WND2B (Wood- Associated NAC 
Domain Transcription Factors) and WND6B were func- 
tional orthologues of Arabidopsis SND1 and master 
switches activating secondary wall biosynthesis during 
wood formation in Populus. Recently, XNDl(Xylem 
NAC Domainl) was reported to influence the differentia- 
tion of tracheary elements and xylem development in 
Arabidopsis by negatively regulating terminal secondary 
wall biosynthesis and programmed cell death in xylem 
vessel cells [29]. 

Although several key switches in regulation of second- 
ary wall formation have been found in the model plants 
Arabidopsis and Populus, key regulators in other plants 
and many downstream genes regulated by these tran- 
scriptional factors are still not clear. In order to identify 
key homologue genes and discover the network control- 
ling cell wall biosynthesis, we carried out comparative 
genome analysis of NST, VND and SND genes across 19 
higher plant species. The analysis revealed that the NAC 
domain transcription factors associated with the second- 
ary cell wall evolutionarily first appeared in Selaginella 
moellendorffii. In fact, among the three groups, only 
VND genes were identified in S. moellendorffii, which is 
evolutionarily earlier than the other two groups. Gene 
expression analysis was carried out to analyse the regula- 
tion of NAC genes associated with secondary cell wall 
biosynthesis in different tissues and revealed that several 
of these transcriptional factors were co-regulated. To 
further characterize the candidate genes involved in the 
regulation of secondary cell wall biosynthesis, we per- 
formed a comprehensive co-regulatory network analysis 
and discovered that some secondary wall-associated 
NAC genes and MYB genes were co-regulated. In 
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addition, co-regulatory network analysis also revealed 
many novel genes and pathways that may be involved in 
cell wall biosynthesis and regulation. 

Methods 

Sequence retrieval and phylogenetic analysis 

Protein sequences and DNA binding domain alignment 
of the NAC transcriptional factor gene family were 
downloaded from Plant Transcriptional Factor Database 
http://planttfdb.cbi.pku.edu.cn. Multiple alignments were 
performed using ClustalX (1.83) software, and the 
Neighbour-Joining (NJ) method was used to construct a 
phylogenetic tree. Genes sharing the same clade with 
the NAC genes controlling cell wall composition from 
Arabidopsis were chosen for further study, which 
resulted in 199 proteins across 19 species. 

Microarray data analysis 

The expression profiling data was acquired from local 
and publically available databases (e.g. GEO and AtGen- 
Express). The signal intensity for each probe set of each 
GeneChip was extracted by Affymetrix software GCOS 
(MAS 5.0). 

Eisen's cluster software http://rana.lbl.gov/EisenSoft- 
ware.htm was applied for cluster analysis. The signal 
intensities of microarray experiments were directly used 
for hierarchical clustering analysis. We employed stan- 
dard data adjustment and SOM (Self-Organizing Map) 
clustering in precedence of hierarchical clustering to 
achieve a better grouping result. 

Gene ontology (GO) analysis was performed for differ- 
entially expressed genes using the EasyGO web server 
[30]. During GO processing, the statistical test method 
used was the chi-square test with FDR p-value < 0.05 as 
the cut-off. 

The gene network data was constructed using Pathway 
Studio http://www.ariadnegenomics.com/products/path- 
way-studio/, ATTED http://atted.jp/ and Hans J. Boh- 
nert's paper [31], and the map was constructed using 
Pathway Studio (version 6.2). 

Plant materials 

Seven tissue samples of rice {Oryza sativa subsp. japo- 
nica var. Nipponbare) were selected for real-time RT- 
PCR (reverse transcription polymerase chain reaction) 
analysis. Rice calli were cultured in N6 solid medium 
[32] and harvested after one month of induction. Root 
samples were harvested from rice seedlings that were 
cultured in a growth container for two weeks. The other 
five samples (penultimate leaf, flag leaf, spikelet, seed 
and stem) were harvested from rice plant grown for 
about four months under natural conditions in Beijing, 
China. 



RNA isolation and real-time RT-PCR 

All rice tissue samples were homogenized in liquid 
nitrogen before isolation of RNA. Total RNA was iso- 
lated using TRIZOL reagent (Invitrogen, CA, USA) and 
purified using Qiagen RNeasy columns (Qiagen, Hilden, 
Germany). Reverse transcription was performed using 
Moloney murine leukemia virus (M-MLV; Invitrogen). 
The cDNA samples were diluted to 8 ng/uL for real- 
time RT-PCR analysis. 

For real-time RT-PCR, triplicate assays were per- 
formed on 1 uL of each cDNA dilution using the SYBR 
Green Master Mix (Applied Biosystems, USA, PN 
4309155) with an ABI 7900 sequence detection system 
according to the manufacture's protocol (Applied Bio- 
systems). The gene-specific primers were designed by 
using PRIMER3 http://frodo.wi.mit.edu/primer3/input. 
htm. The amplification of 18S rRNA was used as an 
internal control to normalize all data (forward primer, 
5'-CGGCTACCACATCCAAGGAA-3' and reverse pri- 
mer, 5'- TGTCACTACCTCCCCGTGTCA-3'). Gene- 
specific primers are listed in Supplemental Table 1 in 
additional file 1. The relative quantification method 
(AACT) [33] was used to evaluate quantitative variation 
between replicates examined. 

Results and discussion 

Identification of genes of NAC transcriptional factors 
controlling the cell wall composition across different 
species 

Plant NAC domain proteins are one of the largest group 
of plant-specific transcription factors. This study mainly 
focuses on SND, NST and VNDs. Phylogeny analysis was 
used to identify the SND, NST and VND members in dif- 
ferent species. The NAM domains of all NAC proteins 
were extracted from the Plant Transcriptional Factor 
Database http://planttfdb.cbi.pku.edu.cn. Multiple align- 
ments and phylogenetic tree analysis showed that the 
genes shared the same clades with Arabidopsis SND, NST 
and VND proteins and thus were used for further analysis. 
This yielded a total of 199 NAC proteins across 19 species 
that possibly regulate secondary cell wall biosynthesis. The 
numbers of potential SND, NST and VND proteins with 
the evolutionary relationship of different species are sum- 
marized in Figure 1, while the number of genes in different 
species may be limited by the Plant Transcriptional Factor 
Database or sequencing information. For example, 
MtNSTl was not included in the Plant Transcriptional 
Factor Database and so was omitted in this study. 

The NAC gene family associated with the secondary 
cell wall biosynthesis evolutionarily first appeared in S. 
moellendorffii (Figure 1). Among the three groups, only 
VND proteins appeared in S. moellendorffii, which was 
evolutionarily earlier than the other two groups. 
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Table 1 Gene Ontology enrichment analysis for the gene list in SNDs-related network 
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Figure 1 Phylogenetic tree of 19 different species and summary of number of SND, NST and VND proteins controlling the cell wall 
composition across these species. 



Evolutionary relatedness of SND, NST, VND genes in 
different species 

A phylogenetic analysis of the 199 NAC proteins was 
performed with ClustalX (1.83) using multiple alignment 
results. A radial NJ tree was generated (Figure 2), and the 
phylogeny tree was divided into three classes differen- 
tiated by specific conserved domains and highlighted 
with three different colours. The three subgroups 
(described in Figure 2) were named NST, SND and VND. 
The orthologues of NSTs (NST1, NST2 and NST3) 
included 48 proteins from 17 species. Through phylogeny 
analysis, NST3 (also named as SND1) orthologues were 
evolutionally closed to NST subfamily. The SND gene 
subfamily encompasses 45 proteins from 17 species. This 
subfamily can also be divided into two subunits: the 
SND2 and SND3 homologues. In the SND3 homologue 
subunit, all genes were from dicots; whereas the SND2 
homologues were from both dicots and monocots. Com- 
pared to SND2 and SND3, the VND subfamily was much 
larger and included four subunits. The first subunit 
included homologues of VND7, which were all from 
dicots. Remarkably, the second subunit only contained 
five proteins from S. moellendorffii. In the other two 



subunits, VND orthologues demonstrated an interspersed 
distribution from different species. 

Expression patterns of SND, NST and VND genes in 
Arabidopsis 

To better understand the functions of SND, NST and 
VND genes, we analysed the gene expression profiles of 
the 12 NAC genes from Arabidopsis using microarray data, 
which included eight diverse Arabidopsis tissues with tripli- 
cate samples using the Arabidopsis ATH1 array [34]. Hier- 
archical cluster analysis was conducted and genes with 
similar expression patterns were clustered together. The 
five NST and SND genes showed a highly similar expres- 
sion pattern in the tissue specific experiments (Figure 3). 
They were preferentially expressed in stems, especially the 
second internodes, and were relatively less expressed in 
leaves, flowers, roots, hypocotyls and mature pollens. In 
siliques, these genes were expressed at an average level. In 
contrast, the VND genes exhibited very different expres- 
sion profiles in different tissues, except for VND6 and 
VND7 which were preferentially expressed in the stem. For 
example, VND1 showed expression in flowers, roots, hypo- 
cotyls and first nodes of the stem while VND3 and VND4 
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Figure 2 Phylogenetic tree of SND, NST and VND proteins across the 19 different species The three clades are indicated in specific 
colours. Members of SND, NST, and VND protein from Arabidopsis are denoted in red. The detailed information of these proteins was listed in 
Supplemental Table 3 in additional file 3. 



were mainly expressed in mature pollens, and second inter- 
nodes of the stem. In general, all seven genes were 
expressed in hypocotyls. With the exception of VND3 and 
VND4, the other VND genes were highly expressed in 
roots. VND2 was also highly expressed in mature pollen 
together with VND3 and VND4. Although the three classes 
of NAC genes showed different expression profiles, they all 
showed lower expression levels in leaves. These microarray 
results suggest that VND genes may have diverse functions 
in different developmental stages or cellular processes. The 



cluster analysis also revealed that several of these transcrip- 
tional factors were co-regulated. 

Expression patterns of SND, NST and VND orthologue 
genes in rice 

Rice is a model monocot plant for molecular and genetic 
studies. Recently, large sets of rice tissue-specific micro- 
array data have been available in the GEO database. 
GSE19024 [35] was selected, including 39 tissues cover- 
ing the entire tissue culture process and life cycle, from 
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Figure 3 Cluster analysis of SND, NST and VND genes of Arabidopsis The gene names are indicated on the right; tissue types are indicated 
above each column. High, average and low levels of expression in a specific tissue are indicated by red, white and blue, respectively. 



rice indica var. Minghui 63, and the gene expression pat- 
terns analysed for rice SND, NST and VND orthologue 
genes that possibly control cell wall composition. Using 
an approach similar to that employed for Arabidopsis 
microarray analysis, four rice NST and SND orthologue 
genes showed similar expression pattern in the tissue 
specific profiles (Figure 4). These four genes were signifi- 
cantly highly expressed in stem tissues, especially in the 
stem_2 (rice stem during heading stage) and were rela- 
tively less expressed in the leaves, endosperm, plumule, 
panicle and callus. The rice VND orthologue genes 
exhibited much more diverse expression patterns in 
these rice tissues. For example, LOC_Os06g04090 
showed preferential expression in stem_2, callus and 
radicle. LOC_Osl0g38834 was highly expressed in pani- 
cle, callus and radicle. We also analysed the expression 
profiles of rice SND, NST and VND orthologue genes in 
the developmental transcriptomes of var. Zhenshan 97, 
another rice indica variety, which showed similar results 
(data not shown) to those of var. Minghui 63. Real-time 



RT-PCR experiments were used to further validate the 
microarray data for those rice NAC genes. Six SND, NST 
and VND orthologue genes were selected to examine the 
relative expression intensity in seven rice tissues, includ- 
ing callus, root, penultimate leaf, flag leaf, spikelet, seed 
and stem. The real-time RT-PCR results (Figure 5) 
mostly matched the microarray results. The expression 
pattern analysis of rice SND, NST and VND orthologue 
genes suggests the existence of conservation of tissue 
specificity in transcription levels between rice and Arabi- 
dopsis, especially for SND and NST orthologue genes 
with preferential expression in developing stem tissues. 

Co-regulatory network analysis for secondary cell wall 
biosynthesis NAC transcriptional factors SND1, SND2 
and SND3 

The Arabidopsis tissue experiment clustering results 
showed that the SND genes were closely clustered 
together and displayed a specific expression in second 
internodes of the stem. The co-expression of genes 
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Figure 4 Cluster analysis of rice SND, NST and VND orthologue genes. The gene names are indicated on the right; tissue types are 
indicated above each column. High, average and low levels of expression in a specific tissue are indicated by red, white and blue, respectively. 



indicates that they may be co-regulated and involved in 
similar molecular regulatory pathways. We applied 
ATTED-II http://atted.jp/ and text mining results to 
predict the possible network where these three NAC tran- 
scriptional factors are involved in secondary cell wall bio- 
synthesis. ATTED-II was built based on 1388 GeneChip 
data in Ambidopsis. The conditions were very diverse, 
including multiple tissue types, abiotic and biotic stress, 
hormone treatments and gene mutants [36]. Using Path- 
way Studio, we built a co-expression gene network for the 
three SND genes. Each SND gene was co-expressed with 
multiple other genes (Figure 6). The three SND genes co- 
expressed with some transcription factor genes, especially 
the MYB genes, e.g. MYB46, MYB52, MYB54, MYB58, 
MYB63, MYB69, MYB85 and MYB 103. These MYB tran- 
scriptional factors are important in modulating cell wall 
biosynthesis [22]. For example, MYB46 is a direct target of 
SND1 and a key player in the transcriptional network 
involved in the regulation of secondary wall biosynthesis 



in Arabidopsis [37]. MYB58 was reported to directly acti- 
vate the expression of lignin biosynthetic genes and a sec- 
ondary wall-associated laccase (LAC4) gene. MYB58 and 
MYB63 are transcriptional activators of the lignin biosyn- 
thetic pathway during secondary cell wall formation in 
Ambidopsis and their expression was regulated by SND1 
[38]. The Class II KNOX gene KNAT7 was also co- 
expressed with SND genes. KNAT7 was shown to nega- 
tively regulate secondary wall formation in Arabidopsis 
and is functionally conserved in Populus [39]. Some other 
transcription factors were also involved in this network, 
e.g. zinc finger, ERF38 and WRKY. 

Furthermore, some secondary cell wall metabolism- 
related genes were co-expressed with SND genes, such as 
LAC genes (LAC2, LACS, LAC10, LAC12 and LAC17), IRX 
(IRREGULAR XYLEM) genes (IRX1, IRX3, IRX6, IRX9, 
IRX12 and IRX 14), CESA4 (CELLULOSE SYNTHASE A4) 
and pectinase related protein. The knockout mutant of 
LAC2 had been reported to reduce root elongation under 
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Figure 5 Real-time RT-PCR analysis of rice SND, NST and VND orthologue genes in different tissues. NAC gene names are indicated on 
each chart; tissue names are indicated on the x-axes. The error bars represent the standard deviations of three replicates. 



PEG-induced dehydration [40] and LAC1 7 mutants appear 
to have a reduced deposition of G lignin units [41]. IRX1 
encodes a member of the cellulose synthase family [42-44], 
IRX3 encodes a xylem-specific cellulose synthase [45], 
IRX6 encodes a member of the COBRA family (similar to 
phytochelatin synthetase) [46], IRX9 encodes a putative 
family 43 glycosyl transferase [47,48], IRX12 (also known 
as LAC4) appears to have laccase activity [41], IRX14 
encodes a putative family 43 glycosyl transferase [48,49], 
CESA4 encodes a cellulose synthase [46,50], and all these 
genes are involved in secondary cell wall biosynthesis. 

Interestingly, several RIC (ROP-INTERACTIVE CRIB 
MOTIF-CONTAINING PROTEIN) genes were also co- 
regulated with SND genes, e.g. RIC2 (involved in pollen 
tube growth and function [51]) and RIC4 (interacts with 
ROP2 during pavement cell morphogenesis and with 
ROP1 to promote apical F-actin assembly [52]). In addi- 
tion, the co-regulatory network analysis also revealed that 
many novel genes were co-expressed with SNDs. 

There were a total of 134 genes involved in this net- 
work (Supplemental Table 2 in additional file 2). GO 
analysis [30] was also performed for these 134 SND co- 
regulated genes (Table 1) to decipher the possible biolo- 
gical pathways in which these genes were involved. Of 
the 134 genes queried, there were 131 genes with anno- 
tated GO items. We used 0.05 as the cut-off of FDR 



adjusted p-value. The most significantly enriched GO 
terms were 'secondary cell wall biogenesis process' 
(GO:0009834, FDR p-value = 2.40E-22), 'cellular poly- 
saccharide biosynthetic process' (GO:0033692, FDR 
p-value = 4.30E-07), 'phenylpropanoid biosynthetic 
process' (GO:0009699, FDR p-value = 2.50E-06), 'tran- 
scription factor activity' (GO:0003700, FDR p-value = 
3.00E-49) and 'laccase activity' (GO:0008471, FDR 
p-value = 3.00E-09). The GO terms related to other bio- 
logical processes were also enriched, e.g. 'response to 
stress' (GO:0006950, FDR p-value = 5.60E-03), 'oxidore- 
ductase activity' (GO:0016491, FDR p-value = 4.00E-04), 
'transition metal ion binding' (GO:0046914, FDR p-value = 
3.50E-05) and transporter activity. Also, most genes were 
localized in nuclear and membrane parts. 

Co-regulatory network analysis of SNDs and GO 
enrichment analysis indicated that most co-expressed 
genes were involved in secondary cell wall biogenesis, 
while we also found that some oxidoreductase activity 
and phenylpropanoid biosynthesis pathway genes were 
co-expressed with SND genes, e.g. peroxidase 64, 
NADPH oxidase and FLS2. Some processes such as car- 
bohydrate synthesis and transport were coordinately 
regulated toward cell wall biosynthesis. There may be 
cross-talk between secondary wall biosynthesis and 
other biological processes. 
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Conclusions 

Combining the bioinformatics data mining and bench 
validation approach, we analysed the NST, VND and 
SND genes across plant species. The comparative geno- 
mic analysis revealed that the group VND of the NAC 
gene family evolutionarily first appeared in S. moellen- 
dorffii. The Arabidopsis and rice gene expression analy- 
sis showed the specific patterns of these NAC genes and 
the conservation of SNDs and NSTs in Arabidopsis and 
rice, and they were preferentially expressed in stems. 
The gene network analysis of SND genes in Arabidopsis 
showed that three SND genes (SND1, SND2 and SND3) 
co-expressed with multiple transcription factor genes, 
especially MYB genes and KNAT7, which are important 



in modulating cell wall biosynthesis. Additionally, the 
co-regulatory network analysis revealed many novel 
genes and pathways that could potentially be involved in 
cell wall biosynthesis and regulation. Nevertheless, there 
may be cross-talk between secondary wall biosynthesis 
and other biological process, such as stress response. 

In summary, these results provided new insight into 
the evolution and the gene regulatory network of a sub- 
group of the NAC gene family controlling cell wall com- 
position from the perspective of bioinformatics. These 
may help us to better understand the possible molecular 
mechanism underlying the regulation network of sec- 
ondary cell wall biosynthesis and, therefore, improve 
plant biomass production. 
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